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Abstract: An in-depth knowledge of the microbiota and metabolites in the solid-state fermentation 
(SSF) of Post-fermented Pu-erh tea (Pu-erh Shucha, PFPT), a Chinese traditional tea with various 
health benefits, is essential to develop modern fermentation technology. In this work, the microbial 
diversity and succession in two laboratory-developed SSF protocols for PFPT were investigated 
using pyrosequencing analyses of the bacterial 16S rRNA and fungal 18S rRNA genes. The active 
bacteria in the initial stages of SSF (seven days) were from the raw materials and environment, with a 
dominance of Proteobacteria in both the raw materials and SSF after seven days. The environmental 
bacteria were inoculated into the tea mass throughout the fermentation process and multiplied, with 
a dominance of Firmicutes at day 14 and 21, and then Firmicutes and Actinobacteria at the last stages 
of fermentation (day 28 and 35). The dominant fungi came from the raw material and were identified 
at the genus level as Aspergillus throughout the SSF process. The contents of tea polyphenols, free 
amino acids, gallic acid, theaflavin, thearubigin, and catechins decreased significantly (p < 0.05), while 
the level of theabrownin increased significantly (p < 0.05). The caffeine content showed no significant 
change (p > 0.05). In total, 30 bacterial and three fungal genera showed significant correlations to 
1-8 and 34 identified tea compounds, respectively (p < 0.05). The dynamics of the microbiota and 
chemical compounds, and correlations between their changes in the SSF of PFPT were revealed, and 
present a foundation for further studies on the microbial effects on chemical compounds. 
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1. Introduction 


Post-fermented Pu-erh tea (Pu-erh Shucha, PFPT) is a well-known traditional Chinese tea with 
special sensory characteristics including mellow taste, stable flavor, and brownish-red color [1]. 
Additional, in vivo, in vitro, and clinical studies have suggested that PFPT has multiple health benefits, 
such as hypolipidemic, antiobesity, antimutagenic, antioxidative, antitumor, free radical scavenging, 
and toxicity-suppressing activity [2]. As a healthy beverage, it is popular in Southeast Asia and has 
also been introduced to the Western marketplace in recent years. 

Unlike the raw material of black tea, fresh tea leaves [3,4], the raw materials of PFPT are sun-dried 
green tea (a variety of green tea products) produced from Camellia sinensis var. assamica (JW Masters) 
Kitamura [5]. In PFPT manufacturing, the post-fermentation process is crucial. In this process, the 
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sun-dried tea leaves are stacked for a few weeks, leading to a series of oxidation, condensation, and 
degradation reactions in the chemical constituents of the tea [6]. From the microbiological view, 
the post-fermentation of PFPT is a solid-state fermentation (SSF). Unlike the chemical changes of black 
tea in fermentation, which are catalyzed by endogenous enzymes in tea leaves [3,4], the chemical 
changes in the SSF of PFPT are catalyzed by the microbial enzymes [1]. Presently, the SSF of PFPT is a 
natural fermentation without the use of starter cultures or the sterilization of raw materials. This gives 
rise to the growth of diverse microorganisms that may be primarily involved in the development 
of the different tastes, flavors, and health benefits of PFPT. Therefore an in-depth knowledge of the 
diversity and dynamics of the microbiota and changes in tea compounds in the SSF of PFPT is essential 
to understanding the formation of the characteristic properties of PFPT, and to develop industrial 
starters and control systems. 

The diversity of microbiota in SSF was investigated by culture methods and many microbial 
species belonging to a variety of families have been identified, such as Aspergillaceae, Bacillaceae, 
Chaetomiaceae, etc. [7]. Additionally, a high level of microbiota diversity was gained using 
culture-independent approaches, such as denaturing gradient gel electrophoresis (PCR-DGGE) [8-11], 
a 16S rRNA gene clonal library [7], and pyrosequencing analyses [12,13]. Recently, we investigated the 
microbial communities and enzymes of an SSF processing sample of PFPT (collected on day 21) using 
454 pyrosequencing and LC-MS/MS approaches [13]. The microbiome and 25 toxic metabolites in 
fermented Pu-erh tea were investigated by Illumina MiSeq platform (Illumina, Inc., SanDiego, CA, 
USA) and a quantitative multiplex metabolite analysis [14]. However, the dynamics of the microbiota 
in SSF had only investigated by PCR-DGGE [8,10,11] and a 165 rRNA gene clonal library [7]. Thus, 
further studies using Next-Generation Sequencing were necessary. Additionally, several reports, 
such as those by Gong [15], Lee [16], and Qin [17], evaluated the changes in the chemical compounds 
during the SSF of PFPT. None of these studies established the influence of the successive microbial 
populations on the changes in the tea’s chemical components during the SSF of PFPT. 

In this work, the microbial succession and dynamics of the tea’s characteristic chemical 
compounds in the SSF of PFPT were studied, and the relationship between microbiota and changes in 
tea compounds was surveyed. This research provides a foundation for further comprehensive studies 
on the formation of the special properties of PFPT and for the development of industrial starters. 


2. Materials and Methods 


2.1. Pu-erh Tea Fermentation and Sample Collection 


The SSF of PFPT in this work was developed according to the traditional method, which is a 
natural fermentation. The raw materials, water, utensils, and environment were not sterilized, and a 
starter was not used. Sun-dried green tea, used as the raw material for the SSF of PFPT, was purchased 
from Puer City, Yunnan Province, China. Two repeated SSFs were performed in our laboratory at 
Kunming City, Yunnan Province, China. A 30-kg sample of the sun-dried green tea leaves was mixed 
with 15 L of tap water to produce a solid content of ~65% (w/v).The leaves were mixed every seven 
days to ensure homogeneity, and tap water was added to keep the solids content at 65%-—75% (as judged 
by the experience of the manufacturer). When the fermented tea mass was reddish-brown and free 
from an astringent taste (~35 days), the tea leaves were air-dried to a water content of ~10%, and 
the fermentation process was stopped. Samples of tea leaves were collected every seven days and 
divided into two parts: one was air-dried and subjected to sensory evaluation according to the protocol 
described by GB/T 23776-2009 [18] and analyses of the tea’s chemical compounds, and the other was 
stored at —80 °C for 454 pyrosequencing analyses. 


2.2. Analysis of PFPT’s Chemical Compounds 


The contents of the tea polyphenols (TPs) and free amino acids (FAA) were determined by 
the spectrophotometric method based on FeSO, and the ninhydrin assay method described by 
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Liang et al. [19], respectively. The main tea pigments, including theabrownin (TB), theaflavin 
(TF), and thearubigin (TR), were analyzed using the spectrophotometric method described by 
Wang et al. [20]. The contents of gallic acid (GA), caffeine (CAF), 1,4,6-tri-O-galloyl-B-D-glucose 
(GG), and catechins, including (+)-catechin (C), (—)-epicatechin (EC), (—)-epigallocatechin (EGC), 
(—)-epicatechin 3-O-gallate (ECG), and (—)-epigallocatechin 3-O-gallate (EGCG) in tea leaves was 
determined by high-performance liquid chromatography (HPLC) using an Agilent 1200 series HPLC 
system consisting of an LC-20AB solvent delivery unit, a SIL-20A autosampler, a CTO-20A column 
oven, a G1314B VWD (280 nm), and an LC Ver1.23 workstation (Agilent Technologies, Palo Alto, CA, 
USA). Separation was completed using the TSK-GEL ODS-80TM column (4.6 mmi.d. x 250 mm, 
Tosoh, Yamaguchi, Japan). The mobile phases were solvents A (0.05 M H3PO04-H20, 5% CH3CN) and 
B (0.05 M H3PO4-H20, 80% CH3CN). Elution conditions and flow rates were: 0-25 min, solvent A 
was reduced from 95% to 55% and B from 5% to 45% (linear gradient); 26 min, solvent B increased to 
100%; 26-31 min, solvent B was kept at 100%; 31-32 min, solvent A increased from 0% to 95%; and 
from 32-38 min, solvent A was kept at 95%; the flow rate was 1.0 mL/min. The temperature of the 
column oven was maintained at 40 °C. The injection volume was 10 uL. The chemical compounds were 
identified in the tea liquids by comparing the retention times of the peaks with those of the standards. 
Each tea sample was extracted twice and each extraction was analyzed twice. 


2.3. Enumeration of Microorganisms in SSF of PFPT 


Tea leaves (25 g) were suspended in 225 mL of 0.9% NaCl buffer. The liquids were serially diluted, 
and the diluted samples were spread on Rose Bengal Medium Agar (5 g peptone, 10 g glucose, 1 g 
KH2POq, 0.5 g MgSO4:7H20, and 20 g Agar per liter) containing 0.1 g/L chloramphenicol and on 
nutrient agar (NA) (5 g peptone, 3 g beef extract, 5 g NaCl and 20 g agar/L) plates for fungal and 
bacterial enumerations, respectively. The Rose Bengal Medium Agar and NA plates were incubated at 
28 °C for 5 days and 37 °C for two days, and the numbers of microorganisms in fermented tea leaves 
were calculated as colony-forming units (CFU)/g. 


2.4. Pyrosequencing Analysis 


DNA extraction, PCR amplification, amplicon quantitation, pyrosequencing, and data analysis 
were developed at Majorbio Bio-Pharm Technology Co., Ltd. (Shanghai, China) on a Roche Genome 
Sequencer GS-FLX Titanium platform (Roche Diagnostics, Basel, Switzerland) using MOTHUR 
software version 1.34 [21,22]. DNA from fermented tea leaves was extracted using the E.Z.N.A.™ 
Soil DNA Kit (Omega Bio-Tek, Norcross, GA, USA) according to the manufacturer’s instructions. 
The TransGen AP221-02 kit with TransStart Fastpfu DNA polymerase (TransGen Biotech, Beijing, 
China) was used in PCR reactions. To analyze the taxonomic composition of the bacterial and 
fungal community, universal primer pairs of 27F (5'/-AGAGTTTGATCCTGG CTCAG-3’) with 533R 
(5’-TTA CCG CGG CTG GCA C-3’), and ITS1 (5’-TCCGTAGGTGAACCT GCGG-3’) with ITS4 
(5'-TCCTCCGCTTATTGATATGC-3’), incorporating FLX titanium adapters, which targeted the V1-V3 
region of the bacterial 165 rRNA gene, and the ITS region of the fungal 18S rRNA gene, were chosen 
for the amplification and subsequent pyrosequencing of the PCR products, respectively. Each PCR 
reaction consisted of 0.4 wL Fastpfu DNA polymerase, 4 uL five-fold Fastpfu Buffer, 2 uL 2.5 mM 
dNTPs, 0.8 uL of each primer, 10 ng of DNA, and ddH,0O to a final concentration of 20 uL. The 
amplification program for the bacterial 165 rRNA gene consisted of an initial denaturation step at 95 °C 
for 2 min. This was followed by 25 cycles, where a cycle consisted of 95 °C for 30 s (denaturation), 55 °C 
for 30 s (annealing), and 72 °C for 30 s (extension), followed by a final extension of 72 °C for 5 min. The 
amplification program for the fungal 18S rRNA gene consisted of an initial denaturation step at 95 °C 
for 2 min. This was followed by 30 cycles, where a cycle consisted of 95 °C for 30 s (denaturation), 55 °C 
for 30 s (annealing), and 72 °C for 30 s (extension), followed by a final extension of 72 °C for 5 min. 
All PCR products were visualized on agarose gels (2% in TBE buffer) containing ethidium bromide. 
For each sample, three independent PCRs were performed. The triplicate products were pooled and 
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purified using an AxyPrep PCR Clean-up Kit (Axygen Biosciences, Axygen, CA, USA). The DNA 
concentration of each PCR product was determined using a Quant-iT PicoGreen double-stranded DNA 
assay (Invitrogen, Carlsbad, CA, USA) and was quality controlled on an Agilent 2100 bioanalyzer 
(Agilent Technologies, Palo Alto, CA, USA). Amplicon pyrosequencing was performed from the A-end 
using a 454/Roche A sequencing primer kit on a Roche Genome Sequencer GS-FLX Titanium platform 
at Majorbio Bio-Pharm Technology Co., Ltd. 

Data preprocessing was performed mainly using MOTHUR software [21,22]. The chimeric 
sequences were excluded using the chimera.uchime command with default parameters [23]. Sequences 
with similarities >97% were clustered into one operational taxonomic unit (OTU) using MOTHUR. 
The taxonomical assignment of each OTU was performed using the classify.seqs command (Naive 
Bayesian Classifier) against the SILVA 16S rRNA gene database (release 111) at an 80% confidence 
level [24-26]. Community richness and diversity indices (Chaol estimator, abundance-based coverage 
estimator (ACE) and Shannon indices, respectively) and rarefaction curves were obtained using 
MOTHUR [27]. 

To compare microbial community structures based on the OTU composition across all samples, 
a principal component analysis (PCA) was carried out using the SPSS software package (SPSS 19.0, 
Chicago, IL, USA). 

The 454 pyrosequencing data generated for this study were submitted to the Sequence Read 
Archive (SRA) and are available under projects SRR1634067 and SRR1633201. The representative 
165 rRNA gene sequences of isolated bacteria in this study were deposited in GenBank under accession 
numbers KP732727-KP732737. 


2.5. Statistical Analyses 


In the analysis of tea chemical compounds, each tea sample was extracted twice, and each 
extraction was analyzed twice. A one-way ANOVA was used to identify statistical differences, 
followed by the least-significant difference method for paired data. Results are expressed as the 


mean + SD. Pearson correlation coefficients between tea chemical compounds and abundant genera 
were calculated. All of the data were analyzed using SPSS 19.0 software packages (SPSS Inc., Chicago, 
IL, USA), and p < 0.05 was considered to be significant. 


3. Results and Discussion 


3.1. SSF of PFPT, Sample Collection, and Variations in Chemical Compounds 


To better understand the natural SSF of traditional Chinese PFPT, laboratory fermentation was 
developed and repeated twice. Every seven days, tea leaf samples were collected and subjected 
to sensory evaluation (Figure 1). After SSF, the tea leaves became dark, the tea liquid became 
reddish-brownish, and the fermented tea leaves became softer than those of the raw material. 
The contents of chemical compounds of tea leaves are listed in Table 1. After SSF, the levels of 
TPs, FAA, EGC, EC, EGCG, GG, ECG, GA, TF, and TR decreased significantly (p < 0.05), and contents 
of TB increased significantly (p < 0.05). These changes in the chemical compounds during the SSF of 
PFPT were in accord with the previous reports of Gong [15], Lee [16], and Qin [17]. It was hypothesized 
that microbiota led to these chemical changes. However, the influence of microbial populations on 
the changes of tea chemical components during SSF of PFPT is not established. So we developed the 
further research to survey the diversity and dynamic of microbiota and change of tea compounds. 
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Table 1. Variation of levels of chemical compounds during solid-state fermentation (SSF) of Post-fermented Pu-erh tea (Pu-erh Shucha, PFPT). 


Sample of Tea Leaves Collected During SSF of PFPT 


Chemical Te 
Variation Trend 
Compounds 0 Day 7 Day 14 Day 21 Day 28 Day 35 Day 
TPs (%) 37.8641.72% 31554264 22634147° 15.7140.674 12 + 0.33 ¢ 10.56 + 0.78 © decreasing (p < 0.05) 
FAA (%) 2.52 + 0.024 2+ 0.03» 1.47 + 0.15° 1.01 + 0.09 4 0.93 +0.15¢ 0.62 + 0.23 ° decreasing (p < 0.05) 
EGC (mg/g) 29.60 £7.94 47.41+8.28>> 33.6546.67% 15.8946.06° 11464251 6974+1.744 — increasing (0-7 day, p < 0.05), decreasing (7-35 day, p < 0.05) 
C (mg/g) 9.83 + 0.763 5.68 +1.17> 4.49 + 0.42 ¢ 3.02 + 0.774 2.02 + 0.71° 1.69 + 0.60° decreasing (0-21 day, p < 0.05) 
EC (mg/g) 13.704 0.12 27.48 +0.62> 17.54 + 4.36 © 11.44 + 2.224 6.05 + 0.944 2.84 + 0.78 © increasing (0-7 day, p < 0.05), decreasing (7-35 day, p < 0.05) 
EGCG (mg/g) 59.73 41.82% 115443575 71542.98¢ 2.86 + 0.93 4 1.39 + 0.614 1.04 + 1.414 decreasing (0-21 day, p < 0.05) 
GG (mg/g) 4.19 + 0.22 ° 0.95 + 0.34 ° 0.90 + 0.14 0.67 + 0.22 © 0.53 + 0.17 © 0.45 + 0.16 © decreasing (0-21 day, p < 0.05) 
ECG (mg/g) 39.63£135*  6.69+4+2.52> 4.13 + 1.26¢ 2.04+0.83¢  0.72+0324 0.574 0.33° decreasing (0-21 day, p < 0.05) 
GA (mg/g) 3.88 + 0.18? 30.24 + 1.68 19.16 + 1.78 © 8.29 + 1.17 ° 1.96 + 0.35 © 0.25+0.11! increasing (0-7 day, p < 0.05), decreasing (7-35 day, p < 0.05) 
CAF (mg/g) 28.50+0.18%  29.53+40.862 32094652? 30.2543.367 30.00+166%  30.5141.757 no significant change (p > 0.05) 
TF (%) 0.35 + 0.018 0.53 + 0.07» 0.45 + 0.03» 0.34 + 0.08 ¢ 0.25 + 0.08 ¢ 0.22+0.024 increasing (0-7 day, p < 0.05), decreasing (7-35 day, p < 0.05) 
TR (%) 7.01 + 0.084 7.78+0.7° 7.48 + 0.44 5.97 + 0.08 ¢ 5.12 + 0.46 ¢ 2.65 + 0.6° increasing (0-7 day, p < 0.05), decreasing (7-35 day, p < 0.05) 
TB (%) 3.23 + 0.08 4 4.16+ 0.212 7.19+049> 113240876 13.3841.724 15.1442.22¢ increasing (p < 0.05) 


Note: Different lowercase letters in superscript (© 44"4¢) in a same row indicate levels of chemical compounds with statistically significant difference (p < 0.05), by one-way ANOVA. 
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Figure 1. The representative sensory evaluation of tea leaves collected on day 0 (a); 7 (b); 14 (c); 21 (d); 
28 (e) and 35 (f) during solid-state fermentation (SSF) of Post-fermented Pu-erh tea (Pu-erh Shucha, 
PFPT). The sensory evaluation was developed according the Chinese standard (GB/T 23776-2009) as 
follows: 3 g of dried tea leaves were infused in 150 mL boiled water for 5 min; the pre-and post-infusion 
tea leaves and the tea liquid after infusion are shown. 


3.2. Diversity of Microbial Communities during SSF of PEPT 


The average initial bacterial and fungal cell numbers in the raw material were approximately 
5.7 x 10* and 1.3 x 10° CFU/g, which increased to about 5.8 x 10” and 3.2 x 10° CFU/g at seven days 
of fermentation, and then increased to about 10° or 10” CFU/g during SSF (Table 2). 
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Table 2. Summary of pyrosequencing analyses of the bacterial 165 rRNA and fungal 18S rRNA genes. 
Sample Cell Number Neve ane: 0 euceee No. of Taxonomic Groups Identified Richness Estimators DIVERSITY Indices 
Reads OTUs (%) 
Means of CFU/g Phylum Family Genus ACE Chao1 Shannon-Weaver Index Simpson 
Bacteria 
D-0-1 : 5409 1338 81.94 7 44 80 8429 4495 4.84 0.07 
D-0-2 5.7 x 10 5529 1482 78.60 8 55 106 9985 5343 5.30 0.04 
D-7-1 5 4637 1479 75.57 14 90 158 10,208 5487 552 0.03 
D-7-2 58 x 10 6691 1988 77.18 10 84 11 15,638 7076 5.50 0.04 
D-14-1 : 4539 1099 84.27 7 60 95 4567 2999 5.43 0.03 
D-14-2 17 x 10 5393 2307 69.63 rp) 86 152 12,646 6382 6.91 0.01 
D-21-1 ; 3472 1302 73.30 5 38 63 7453 3622 5.90 0.03 
D-21-2 12 x 10 5811 2071 75.17 10 59 92 10,983 6042 6.57 0.01 
D-28-1 P 3175 1261 71.28 j 53 76 7165 3621 5.97 0.02 
D-28-2 2.0 x 10 3649 1178 76.38 7 57 85 6782 3989 5.34 0.04 
D-35-1 ; 3119 1138 73.45 5 51 81 5953 3319 5.66 0.02 
D-35-2 8.3 x 10 3137 948 78.39 7 52 77 5157 2687 4.96 0.06 
Fungi 

D-0-1 . 2189 162 96.62 5 41 58 330 258 a 0.26 
D-0-2 13 x 10 1395 92 97.56 5 27 36 127 127 227 0.31 
D-7-1 : 6305 32 99.71 3 3 4 128 63 0.4 0.84 
D-7-2 3.2 x 10 5773 30 99.65 Z 3 4 261 125 0.29 0.89 
D-14-1 2 5718 41 99.56 2 3 3 84 116 0.98 0.48 
D-14-2 9.0 x 10 5737 33 99.65 e 3 3 293 128 0.65 0.71 
D-21-1 : 6211 41 99.60 2 2 3 270 101 0.50 0.79 
D-21-2 9.6 x 10 5487 61 99,42 2 3 4 93 123 1.10 0.47 
D-28-1 : 5314 66 99.36 2 re 13 80 136 0.66 0.78 
D-28-2 L210 5955 49 99.63 4 3 4 114 67 0.66 0.75 
D-35-1 F 8050 99 99.30 2 6 9 325 195 0.78 0.75 
D-35-2 6.6 x 10 5199 67 99,29 3 4 5 248 123 0.87 0.69 
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The microbial communities were analyzed using a pyrosequencing-based analysis of the bacterial 
16S rRNA gene and fungal 185 rRNA gene in the 12 samples of fermented tea leaves. After removing 
low-quality and chimeric sequences, a total of 54,561 bacterial and 63,333 fungal high quality sequences 
with average read lengths of 472 and 497 bp, respectively, were obtained. Each bacterial library 
contained 3119-6691 reads, with different phylogenetic OTUs ranging from 948 to 2307, which were 
classified into 5-14 phyla, 38-90 families and 63-158 genera (Table 2). In them, 15,490 sequences 
were classified as cyanobacteria, which are probably chloroplast sequences that originated from the 
tea plant. These sequences were removed prior to further analysis. Each fungal library contained 
1395-8050 reads, with different phylogenetic OTUs ranging from 30 to 162, and was classified into 
2-5 phyla, 2-41 families, and 3-58 genera (Table 2). 

Rarefaction curves at a 3% cutoff are shown in Figure 2, and the Chaol estimator, abundance-based 
coverage estimator (ACE), and Shannon index are listed in Table 2. The bacterial diversity indices 
were significantly greater than that of fungi (p < 0.05), the bacterial rarefaction curve did not approach 
a plateau at a similar level, and more bacterial taxonomic groups were identified, demonstrating that 
bacterial communities in fermented tea leaves have a high diversity level. The fungal diversity indices 
in raw material (day 0) was significantly greater than that of fermented tea leaves (day 7-35) (p < 0.05), 
indicating that a decrease of fungal diversity during SSF. This decrease of fungal diversity was also 


reported by Zhang et al., recently [14]. It seems that as soon as water is added to start the fermentation, 
the fungal diversity decreased. 
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Figure 2. Rarefaction curve analyses of the 165 rRNA gene (a) and 185 rRNA gene (b) of the bacterial 


and fungal communities, respectively, present in the SSF of PFPT. Operational taxonomic units (OTUs) 
are defined as sequences sharing 97% nucleotide sequence similarity. 
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3.3. Dynamics and Source of Bacteria during SSF of PFPT 


A heat map of major bacteria (>1%) in tea leaves at the genus level is given in Figure 3a. A PCA 
analysis is shown in Figure 3b. The bacterial community in the raw material (day 0) was different from 
those in other samples, with dominant bacteria of Sphingobacterium (16.48%) and Methylobacterium 
(16.39%). After seven days of fermentation, the dominant bacteria had changed to Pantoea (21.32%) and 
Pseudomonas (18.87%). In the 14 days of SSF, the genera of Paenibacillus and Bacillus were dominant and 
accounted for 26.21% and 16.73% of total reads, respectively. Achromobacter (17.62%) was dominant 
at 21 days. In the sample collected on day 28 of SSF, the dominant bacteria were identified as 
Staphylococcus (20.78%) and no_rank_Dermabacteraceae (16.91%). At the last stage of SSF (day 35), 
no_rank_Dermabacteraceae (30.82%) were dominant, followed by Staphylococcus (19.62%). At the 
phylum level, the pyrosequencing analysis indicated that the dominant bacteria were Proteobacteria, 
Firmicutes and Actinobacteria. By metagenomic analysis, Lyu showed that dominant bacteria in one 
fermented sample of PFPT were Proteobacteria (23.56%), Actinobacteria (23.35%), and Firmicutes 
(11.37%) [12]. By using a 16S rRNA gene clone library, we also demonstrated that the dominant bacteria 
during SSF of PFPT were Proteobacteria and Firmicutes [7]. Recently, Zhang identified that the most 
commonly observed bacterial taxa in 31 Pu-erh tea samples belonged to Firmicutes, Actinobacteria, 
and Proteobacteria [14]. These samples in different reports have similar dominant bacteria, suggesting 
that Proteobacteria, Firmicutes, and Actinobacteria may be the core bacteria in SSF of PFPT, and they 
may play crucial roles in the quality formation of the PFPT. However, more samples from different 
areas and producers should be investigated to validate this hypothesis. 

To investigate the sources of the bacteria, a series of comparisons of major bacteria (>1%) at the 
genus level were performed (Figure 4). A comparison of raw material and the seven-day fermented 
sample showed that they shared six bacterial genera, which accounted for 67.79% of the total reads of 
seven-day fermentation (Figure 4a). This revealed that partial bacteria in initial stages of SSF came from 
the raw material. Additionally, four genera were newly detected in the seven-day fermented sample. 
This demonstrated that a portion of the bacteria came from the fermentation environment (room, air, 
or water). Comparisons of bacterial genera between the samples of day 7 and 14, day 14 and 21, and 
day 21 and 28 showed that six, six, and five genera, respectively, accounting for 55.49%, 18.61%, and 
26.89%, respectively, of each sample’s reads were newly detected. A comparison of 35-day-fermented 
tea leaves and raw material showed that they shared only 18.35% of bacteria (Figure 4f). Together, these 
comparisons indicated that active bacteria were from the raw material and fermentation environment 
during the initial stages, and the bacteria in the fermentation environment continued to naturally 
inoculate into tea mass during the SSF. This primary comparison suggested that environmental bacteria 
are involved in the fermentation of PFPT; this is possibly the reason for the geographical characteristics 
of PFPT. 
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Figure 3. Heat map of relative abundances of major bacterial genera (>1%) (a); principal component 
analysis (PCA) analysis of bacterial communities in fermented tea leaves during the SSF of PFPT based 
on the relative abundances of the sample reads at the genus level (b). 
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Figure 4. Comparisons of major bacteria (>1%) at the genus level between day 0 and 7 (a); day 7 and 14 
(b); day 14 and 21 (c); day 21 and 28 (d); day 28 and 35 (e); and day 0 and 35 (f) during the SSF of PFPT. 


3.4. Dynamics and Source of Fungi during SSF of PFPT 


In raw material, the dominant fungi were classified as no_rank_Fungi (65.39%). When the water 
content of tea leaves was increased, the fungi belonging to Aspergillus were multiplied to dominate 
throughout the SSF process. They accounted for 94.08%, 74.79%, 79.27%, 93.76%, and 95.79% of the total 
fungal community at day 7, 14, 21, 28, and 35, respectively. Additionally, fungi belonging to Rhizomucor 
multiplied and accounted for 3.73% at day 7, and increased to 24.54% and 20.63% at day 14 and 21, 
respectively, and then decreased to 3.07% at day 35 (Figure 5a). The differences in the dominant fungi 
between the raw material and fermented samples suggested a large shift in the fungal communities, 
which was corroborated by the separation of fungal communities in the raw material (day 0) from 
other fermenting samples (day 7, 14, 21, 28, and 35) in the PCA analysis (Figure 5b). The similarity in 
the dominant fungi in each fermented tea leaf sample suggests that fungal communities did not change 
significantly during the SSF of PFPT, and this was corroborated by the groupings of the fermented 
samples (day 7, 14, 21, 28, and 35) in the PCA analysis (Figure 5b). Fungi belonging to Aspergillus were 
identified as dominant throughout the SSF process of PFPT, and this corroborated previous reports 
that the major fungi involved in the SSF of PFPT were A. niger and B. adeninivorans [8], or A. niger, 
S. cerevisiae, and P. glabrum [11], and several Aspergillus spp. [1]. Comparisons of major fungi (>1%) at 
the genus level showed that they came from the raw material, suggesting that the raw materials were 
responsible for the active fungi. To our knowledge, this is the first work to survey the sources of the 
bacteria and fungi in the SSF of PFPT. 
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Figure 5. Fungal communities in tea leaves at the genus level during the SSF of PFPT (a); PCA 
analysis of fungal communities in fermented tea leaves during the SSF of PFPT based on the relative 
abundances of the sample reads at the genus level (b). Minor groups are composed of the parties 


showing a percentage of reads <1.0% of the total reads in each sample. 


3.5. Correlations among Microbiota and Tea Components 


To explore microbial effects on tea components, correlations between chemical compounds and 
relative abundances of the bacterial and fungal genera were conducted and are shown in Tables 3-5. 
In total, 30 bacterial genera showed significant correlations to 1-8 tea compounds. For example, 
unclassified Actinobacteria have negative significant correlations to the contents of TPs, FAA, EGC, 
C, EC, TF, and TR; and have positive significant correlations to TB (Tables 3 and 4). For fungi, 
Aspergillus showed negative significant correlations to the contents of C, EGCG, ECG, and GG; OTUs 
assigned to no_rank_Fungi have positive significant correlations to the contents of C, EGCG, ECG, 
and GG; and Penicillium had positive significant correlations to the contents of EGCG, ECG, and GG 
(Table 5). The compounds had positive significant correlations to 2-10 bacterial genera and showed 
negative significant correlations to 0-5 bacterial genera (Table 3). This preliminary data suggested that 
microbiota with significant correlations to tea compounds may be involved in the dynamics of chemical 
compounds; however, the roles of these microbiota in the SSF of PFPT need further monoculture 


fermentation studies. 
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Table 3. Significant correlation coefficients between the tea’s chemical compounds and bacterial genera. 
Compounds TPs AA CAF EGC Cc EC EGCG GG ECG GA TF TR TB 
positive Sphi Q : ' a5 B z E : : ‘ “ x 
Sher phingomonas, Sphingomonas, Stenotrophomonas, Erwinia, Sphingomonas, is ax  Sphingobacterium, Sphingobacterium, Sphingobacterium, Erwinia, ae iy as «  Prevotella, 
ae 0.745 * 0.774 * 0.873 * 0.955 * 0.889 * Erwinia, 0.957" 0 ggg 0.901 * 0.901 * 0.948 ** Broa Ooo Eran, 07ES @ O77 * 
Rhizobium, 0.802 * Rhizobium, 0.829 * ee Pantoea, 0.792* — Lactococcus, 0.741*  Pantoea, 0.840 * oe oes ee Aloe Pantoea, 0.765 * coe Kocuria, 0.808 * 
m Herbaspirillum, Enterobacter, ae Enterobacter, Sphingomonas, Sphingomonas, — Sphingomonas, Enterobacter, Streptococcus, Bacteroides, 
Herbaspirillum, 0.774 * 0.797 *P 0.825 * Rhizobium, 0.929 ** 0.912 * laine noose oeates Pantoea, 0.839 * 0.943 * Re * 0.797 * 
Chryseobacterium, Chryseobacterium, Rhodococcus, Herbaspirillum, Rhodococcus, Rhizobium, Rhizobium, Rhizobium, Enterobacter, Rhodococcus, Carnobacterium, pian gue 
0.862 * 0.880 * 0.851 * 0.906 ** 0.848 * 0.992 ** 0.996 ** 0.992 ** 0.967 ** 0.896 ** 0.763 5 923 * oe 
k no_rank_ 
Bes séoceacta Peptostreptoc- Flavobacterium,  Chryseobacterium, Streptococcus, Herbaspirillum, Herbaspirillum,  Herbaspirillum, — Rhodococcus, Streptococcus, 
cage Prococeaceae,  accaceae, 0.826 * 0.966 ** 0.788 * 0.995 ** 0.996 ** 0.997 ** 0.935 ** 0.814 * 
iy 0.821 * 
Kluyvera, Kluyvera, Carnobacterium, oe Flavobacterium, Chryseobacterium, | Chryseobacterium, Chryseobacterium, Flavobacterium, Flavobacterium, 
0.957 * 0.956 * 0.768 * Gene ee 0.902 ** 0.990 * 0.990 * 0.989 * 0.850 * 0.870 * 
Carnobacteriuni no_rank_ no_rank_ no_rank_ no_rank_ no_rank_ no_rank_ no_rank_ 
Lactobacillus, 0.750 * 0.903 ** f Spori-chthyaceae, Kluyvera, 0.970 ** Spori-chthyaceae, Bactero-idales, Bacteroid-ales, Bacteroid-ales, | Sporichth-yaceae, Spo-richthyaceae, 
: 0.866 * 0.891 ** 0.845 * 0.849 * 0.847 * 0.869 * 0.825 * 
no_rank_ no_rank_ 
Carnobacterium, Mesorhizobium,  Novosphingobium, — Mesorhizobium, aoe ankPeplost Peptostre- Peptostre- Mesorhizobium, Mesorhizobium, 
0.919 * 0.894 ** 0.852 * 0.906 ** Dore Ptococcaceae, —_ptococcaceae, 0.942 ** 0.882 * 
, 0.906 ** 0.920 ** 
Lactobacillus, no_rank_ ax  Kluyvera, Kluyvera, 
0.751 * Bacteria, 0.773* — Kluyvera, 0.921 9 aoe ax 0.912 ** 
Carnobacterium, Novosphingobium, Novosphingobium, Novosphingobium, 
0.860 * 0.971 ** 0.984 ** 0.976 ** 
negative 7 Ses no_rank_ 2 7 ae 
ionificant Bacteroides, Bacteroides, Peptost- Prevotella, Bacteroides, Kocuria, 0.777 * Prevotella, Prevotella, Erwinia, 
set fs —0.803 * —0.772* reptococcaceae, —0.753 * —0.763 * Ie —0.771* —0.969 ** —0.740 * 
correlations —0.801 * 
unclassified_ : unclassified_ 5 
Staphylococcus, : : Kocuria, ; : Staphylococcus, Staphylococcus, — Kocuria, Kluyvera, 
i Actinobacteria, .. Act-inobacteria, = z a i" 
—0.749 0.852 * —0.754 0.803 * —0.779 —0.809 —0.957 —0.872 
unclassified_ unclassified_Act- unclassified_Act- . 
; A Staphylococcus, : : : : Staphylococcus, Streptococcus, 
Actino-bacteria, _0.767 * inobacteria, inobacteria, —0.850 * —0.732 * 
0.856 * ; —0.793 * —0.799 * ; : 
unclassified_A- Kocuria, uncl assified_- Lactobacillus, 
ctinobacteria, _0.792 * Actinobacteria, _0.731 * 
—0.815 * : —0.942 ** : 


Carnobacterium, 
—0.908 ** 


Note: * and ** represent significance levels of 0.05 and 0.01, respectively. 
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Table 4. Summary of the significant correlation coefficients between bacterial genus and the tea’s 


chemical compounds. 


Bacterial Genus 


Positive Significant Correlations 


Negative Significant Correlation 


unclassified_Actinobacteria 
Kluyvera 
no_rank_Peptostreptococcaceae 
Carnobacterium 
Erwinia 
Herbaspirillum, Rhizobium, Sphingomonas, Chryseobacterium 
Enterobacter 
Kocuria 
Staphylococcus 
Bacteroides 
Flavobacterium, Mesorhizobium, no_rank_Sporichthyaceae, 
Pantoea, Rhodococcus 
Novosphingobium 
Prevotella 

Streptococcus 
Sphingobacterium 

Lactobacillus 

Microbacterium, no_rank_Bacteroidales 

Lactococcus 

no_rank_Bacteria 
Paenibacillus, Stenotrophomonas 
Pseudomonas 


TB 
TPs, AA, C, EGCG, GG, ECG 
TPs, AA, C, EGCG, GG, ECG 
TPs, AA, EGC, C, TR 
EGC, EC, GA, TF, TR 
TPs, AA, C, EGCG, GG, ECG 
EGC, EC, GA, TF, TR 
TB 


TB 
EGC, EC, GA, TF 


C, EGCG, GG, ECG 
TB 
EC, TF, TR 
EGCG, GG, ECG 
TPs, C 
EGCG, GG, ECG 
c 
EC 
CAF 
GA 


TPs, AA, EGC, C, EC, TF, TR 
TB 


EGC, EC, TF, TR 
TPs, EGC, EC, TF, TR 
TPs, AA,C 


EGC, TF, TR 
TB 


TB 


Table 5. Significant correlation coefficients between contents of the tea chemical compounds and the 


fungal genus. 


Compounds Fungal Genus, Correlation Coefficients 
Cc Aspergillus, —0.875 * No_Rank_fungi, 0.867 * - 
EGCG Aspergillus, —0.959 ** No_Rank_fungi, 0.984 ** Penicillium, 0.863 * 
GG Aspergillus, —0.974 ** No_Rank_fungi, 0.990 ** Penicillium, 0.875 * 
ECG Aspergillus, —0.963 ** No_Rank_fungi, 0.988 ** Penicillium, 0.872 * 


Note: * and ** represent significance levels of 0.05 and 0.01, respectively. 


3.6. Staphylococcus in SSF of Pu-erh Tea 


Staphylococcal food poisoning is a common food-borne disease worldwide following the 


ingestion of staphylococcal enterotoxins that are produced by enterotoxigenic strains, mainly 
Staphylococcus aureus [28]. S. aureus can grow ina wide range of temperatures, pH levels, and sodium 
chloride concentrations, and in many food products [29]. Bacteria belonging to the genus Staphylococcus 
accounted for 20.78% and 19.62% of the total reads at the last stages of SSF (day 28 and 35, respectively), 
suggesting that S. aureus occurs in the PFPT. To verify this possibility, a Baird—Parker plate assay to 
detect S. aureus and the isolation of bacteria in the day 28 and day 35 samples was performed. In total, 
10 and one isolates in which the 165 rRNA genes (KP732727-KP732737) shared 100% similarity to 
those of S. sciuri subsp. sciuri DSM 20345(T) and S. saccharolyticus ATCC 14953(T), respectively, were 
identified by an EzTaxon-e search [30]. However, S. aureus was not detected. The S. sciuri subsp. 
sciuri and S. saccharolyticus present in the fermentation of PFPT need further study to determine 
whether they are harmless symbionts or pathogens. However, their detection only at the last stage 
of SSF by pyrosequencing analysis suggested that the Staphylococcus came from the fermentation 
environment. Thus, we suggest that the safety of the environmental bacteria correlates with the safety 
of PFPT fermentation. 


4. Conclusions 


Spectrophotometry and HPLC revealed that the levels of TPs, FAA, EGC, EC, EGCG, GG, ECG, 
GA, TE, and TR decreased significantly (p < 0.05), while the TB content increased significantly (p < 0.05) 
after SSF. The pyrosequencing analysis and comparisons suggested that the bacteria in the initial stages 
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of the SSF were from the raw material and environment, with Proteobacteria dominant in the raw 
material and at day 7 of SSF. Environmental bacteria were inoculated into the tea mass throughout the 
fermentation process and Firmicutes became the dominant bacteria at day 14 and 21. During the last 
stages of fermentation (day 28 and 35), Firmicutes and Actinobacteria were dominant. The dominant 
fungi were identified as genus Aspergillus in the SSF process. Correlations between microbial succession 
and the dynamics of the chemical compounds were conducted. This work provided information on 
the microbial succession and dynamics of chemical compounds, and the correlations between these 
two changes in the SSF of PFPT. 
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